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Abstract 

Bloch equations give a quantum description of the coupling between an atom and a 
driving electric force. In this article, we address the asymptotics of these equations for high 
frequency electric fields, in a weakly coupled regime. We prove the convergence towards 
rate equations (i.e. linear Boltzmann equations, describing the transitions between energy 
levels of the atom). We give an explicit form for the transition rates. 

This has already been performed in |BFCD03] in the case when the energy levels 
are fixed, and for different classes of electric fields: quasi or almost periodic, KBM, or 
with continuous spectrum. Here, we extend the study to the case when energy levels are 
possibly almost degenerate. However, we need to restrict to quasiperiodic forcings. The 
techniques used stem from manipulations on the density matrix and the averaging theory 
for ordinary differential equations. Possibly perturbed small divisor estimates play a key 
role in the analysis. 
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In the case of a finite number of energy levels, we also precisely analyze the initial 
time- layer in the rate aquation, as well as the long-time convergence towards equilibrium. 
We give hints and counterexamples in the infinite dimensional case. 

Keywords: density matrix, Bloch equations, rate equations, linear Boltzmann equation, 
averaging theory, small divisor estimates, degenerate energy levels. 



1 Introduction 

Bloch equations model the time evolution of a quantum mechanical system described in the 
density matrix formulation and driven by an electromagnetic field. This formalism is very 
precise but sometimes difficult to interpret and to use in practical simulations. It is therefore 
useful to find asymptotic models under appropriate scaling assumptions. In a former article 
[BFCD03 such a program is performed and leads to the rigorous derivation of rate equations, 
which are often used in the physics literature. Here, we want to extend these results to the case 
when the energy levels of the system are almost degenerate. There are many examples of such 
almost degeneracies. This is the case for example of Zeemann hyperfine structures in complex 
molecules, or quantum dots submitted to an external magnetic field. High levels of an atom 
are also almost degenerate since there is an infinite number of levels with accumulation value 
at the ionisation energy. 



1.1 Bloch equations 

According to the quantum theory, matter is described via a density matrix p, whose diagonal 
entry p(t, n, n) is -in the eigenstates basis- the population of the n-th energy level at time t, 
and the off-diagonal entry p(t, n, m) is linked to the transition probability from level n to level 
m (conditioned by the corresponding populations). One may think of a collection of identical, 
uncoupled atoms, with discrete energy levels. We refer the reader to |Boh791 |Boy92| [QTDRG88 , 



ILou91| |NM92| I5SL77| lBid03j for textbooks about wave/matter interaction issues, where Bloch 
equations occur. To treat the mathematical problem, we use a dimensionless version of these 
equations and consider that the density matrix p(t, n, m) is governed by: 

e 2 d t p(t,n,m) = —iuj e (n,m)p(t,n,m)+Q e (p)(n,m) (1) 



V ( — ,n, kj pit, k,m) -V (-^, k,m ) p(t,n, k) 



Integers n G N, m G N, and k G N are labelling discrete energy levels, and time t belongs to 
M + . In the case of a finite number of energy levels, we add the restriction n < N, m < N, 
k < N. 

The dependence of the density matrix on the small parameter e will always be implicit. 
Here, the system is forced by a high frequency electromagnetic wave which is described by its 
oscillatory amplitude 4>{t/e 2 ) and contributes to the quantity 

V ( — ,n,m) = ( -= J V(n,m), 
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where the interaction coefficient V(n, m) 6 C is (up to a rescaling) an entry in the dipole 
moment matrix which is Hermitian: V(m,n) = V(n,m)*. 

The parameter e occurring in the scaling plays two roles. First, the coupling between atoms 
and the wave is small, of order e, and its cumulated effects are considered over long time scales, 
of size 1/e 2 . This is the setting of the so-called weak coupling regime (see |Spo80[ |Spo91[ [vH55 , 
vH57j). Second, the scaling in Eq. (JTJ) produces resonances between the eigenfrequencies u>k/e 
of the free atom and the electromagnetic source <p{t/e 2 ). 

The quantity u £ (n,m) = u(n,m) + e p 5(n,m) is the transition energy between levels n and 
m. In the non-degenerate framework u(n,m) = u(n) — oj(m) 6 1 is the difference between 
the energies u(n) and uj{m) of levels n and m respectively. Here we model almost degenerate 
levels replacing the energy uo{n) by u(n) + e p 5(n) where p > and 5(n) e R. Therefore 
8(n, m) = 5(ri) — 5(m). We notice that a large value for p means that levels n and m are very 
close from one another if u(n, m) = 0. 

Last, relaxation terms are modelled by the operator Q £ which reads 

Q £ (p)(n, m) = — £ M 7(n, m)p(t, n,m), if n ^ m, 

= e 2 ^2 n )p(ti k) - W(n, k)p(t, n, n)]^j , if n — m. 

We assume that the Pauli coefficient W(n, m) is non-negative, as well as the longitudinal coef- 
ficient j(n,m), which is besides symmetric: / -f(n,m) = / ~f(m,n). To simplify further notations, 
we extend its definition to the case when n = m introducing 7(71, n) = 0. 

We assume that < p, < 1/2. This is already the case in the nondegenerate case BFCD03 . 
This threshold value for p arises in the estimates and we do not know whether it is optimal 
or not. However it is no wonder that there is a threshold value. Indeed, in the case when 
p = 0, the initial Bloch equation (0) is time- irreversible and the asymptotic equation that we 
derive in this paper is also time-irreversible. On the other hand, in the opposite case when 
every coefficient 7(71, m) is identically zero (which can also be interpreted as p — 00), the initial 
Bloch equation is time-reversible and the nature of the problem has changed. 

The relaxation operator Q £ determines an equilibrium state, to which the system has a 
tendency to relax when no electromagnetic wave is applied. If the threshold p < 1/2 is not 
necessarily optimal, it is very important that p < 2. This means that the off-diagonal terms of 
the density matrix are more rapidly decreasing that the diagonal terms. We therefore expect 
them to play no significant role on the time scale 1/e 2 . 

Finally, we associate to Eq. (P) an initial datum p(0, n, m) which satisfies 

p(0, n, m) = 0, Vn 7^ m, and p(0, n, n) > 0, Vn. (2) 

The vanishing of the off-diagonal entries means that we are dealing with "well prepared" initial 
data. This is a standard assumption in the field (see e.g. [K L571 ILK581 fZwa66 ). 

Remark. Note that the small coupling in Eq. (0) a priori induces perturbations of size 1/e on 
time scales 1/e 2 . As we shall see, the very Hamiltonian nature of the equation actually makes 
these perturbations of size 1 . 
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1.2 Rate equations 

To describe the asymptotic dynamics, we are only interested in the diagonal part of the density 
matrix (the populations) and set 

Pa(t,n) := p(t,n,n), 

which is the occupation number of the ra-th level. The limiting model we look for as e — > is 
a system of rate equations which are Boltzmann type equations, in the form 

dtPT = (wr d hpT, (3) 

where the approximate populations p°| pp are viewed as a vector, and the modified relaxation 
matrix operator (W / e mod )ti is defined from coefficients (n, k) via the following notation. 

Notation 1. To coefficients A(n, k), n ^ k (which may possibly be time dependent: A(t, n, k) ) 
we associate the matrix A$ (respectively A$(t)) through 

A$(n,k) := A(k,n) ifn^k, 
A$(n, n) := — ^] A(n, m) else. 

If A(n, k) G H l^ll,, A$ is a bounded operator on l p , 1 < p < oo (see Appendix \8 . 1\) . 

Apart from |BFCD03j . in the past few years an extensive attention has been paid on the 



rigourous derivation of Boltzmann type equations from dynamical models of (classical or quan- 
tum) particles or models for the interaction of waves with random media. Convergence results 
in the case of an electron in a periodic box may be found in |Cas99[ ICas02[ ICasOlj . We also 
mention the non-convergence result established in jCP02| IDP03 in a particular, periodic situ- 
ation. For the case when an electron is weakly coupled to random obstacles, the reader may 
refer to |EY0n[ |Spo771 |Spo80[ Spo91| and |KPR 96j for the formal analysis performed. The 



computation of the relevant cross-sections is performed in |Nig.Q.6J . All these results address the 
case of a linear Boltzmann equation. A nonlinear case is studied in [BCEP03 . 

As proved in |BFCD03 , in the non-perturbed case (S(n) = 0) the action of the wave 



transforms the relaxation operator W$ associated with matrix W into a limit relaxation operator 

or-");. 

The goal of this article is to specify the precise form of matrix W £ mod (n, k) for the perturbed 
case for a fixed e. Once this rate equation is derived, we also give results about the limit e — ► 
and the size of the time-layer leading to equilibrium according to the value of p/p. 

1.3 Main results 

1.3.1 The modified relation operator (W r £ mod )ti 
The modified transition rates W r £ mod (n, k) is defined by 

ty £ mod (n, k) := (^ £ ) dom (n, k) + W(n, k), 
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and the non-negative coefficient (\I/ e ) dom (n, k) is given (see Proposition 0] and Theorem EJ) by 

(..>-(„,*, ^2in>a)r-E ^. n)2 + Mn ;^.^ + ^. n)| ^ / if , = 0, 
!=2 ' y t".*)i' ^ (t ,?* , ^ (M) . E W ^>°- 

Here, the electromagnetic wave is determined by its frequency vector c<j £ R r and its Fourier 
coefficients (fi a : 

4>(t) = 4> a exp(ia ■ ut), where we denote a ■ u := a.\Ui + • • • + a r u r . 

Under assumptions that we specify in Section |2J), we prove the following: if p is a solution 
to Eq. (JT|). and p^ pp the solution to Eq. (jSJ) with the same initial datum, then for all T > 0, 
there exists C > such that 

||Pod|U°°([o,T],*i) < Ce 1 '^ and ||p d - Pd PP |U-([o,T],z2) < C(e A< + e 1_2M ). 



1.3.2 The asymptotic state 

We study the convergence of the solution p to Bloch equation (0) on some fixed time interval 
[0,T], as e goes to zero. This dynamics is thus given by the corresponding solution p^ pp to 
Eq. ®. Since the operator (W™ d )§ has a non-positive spectrum (see the Appendix I8.1J1 . one 
may think that, as time grows, p^ pp approaches some equilibrium state, i.e. a stationary state 
p El 2 belonging to the kernel of {W™ A )§. We describe this convergence carefully in the case of 
a finite number of quantum levels. In the case of an infinite number of levels, we exhibit simple 
examples for which convergence does not occur. 



1.3.3 The case of a finite number of levels 

The evolution of pa can be summarized as: 

(i) The levels n which do not resonate with the wave ((\l/ £ ) dom (n, k) = for all k) and are 
not coupled to others via relaxations either (W(n, k) = for all k) remain constant. 

(ii) Polarization: the remaining levels belonging to the space (Z 2 )-orthogonal to the kernel 
of the sing ular part (^ e )^ mg of (^ e ) dom vanish in time 0(6°). The precise definitions of 
(\l/ e )jj ing and a > are given in Sectional 

(iii) After this time, the evolution of the remaining levels is given by an e-independent system 
forlK PP , 

d t u P f p = u(w + ^ onsing )»riPd pp - 

Here, n is the orthogonal projector onto the space orthogonal to the ones treated in 
the first two steps. Also, we set (^ £ ) nonsin g : = (^ e ) dom - (^ e ) sin s, and ^° onsin s is the 
value of (v|/ e } nonsm g as e vanishes. For each initial datum p(0), the solution to this system 
converges, as t goes to infinity, to some (unique) equilibrium state p (see Appendix 18. 1J) . 
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In the unperturbed case 6 = |BFCD03] . we have a = /i. The perturbation can affect the 
exponent a, the transition rates given by (\l/ e )jj mg , as well as the final equilibrium (determined 
by (ty £ )™ nsmg ). The precise value of o and the possibility of modifying the limiting operator 
according to the value of fi/p are given in Tabled 



Li/f) 


a 


difference with the unperturbed case 

|bTCDoh] 


< (i/p < 1 




none 


fjt/p = 1 




transition rates 


1 < /j,/p < 2 


2p-fi 


transition rates and time-layer 


2 < /i/p < oo 




projector IT and asymptotic state 



Table 1: Consequences of the perturbation for a finite number of levels. 



1.3.4 The case of an infinite number of levels 

For an infinite number of levels, Proposition 0] and Theorem |U] hold. They give the approxima- 
tion of the solution to Eq. (jl} by the solution to Eq. 0. The transition rates W(n, k) are again 
transformed into W™ d (n, k). However, the analysis of the asymptotic behavior of the solution 
to Eq. (jHJ) is more intricate than in the finite case. In Appendix 18.1.21 we gi ye examples for 
which there does not exist any possible equilibrium state in I 2 (see Corollary 1 16|) . 

It is still true that non- interacting levels remain constant, as in item (i) of Section 11.3.31 
Nevertheless, we are unable to prove the precise polarization property of the two other items, 
even if we expect that the results of Table Q hold. An equivalent to Table Q for an infinite 
number of levels is given in Section IFTHj for the values of /i/p for which we can conclude. 



1.4 Outline 

The article is organized as follows. Section El is devoted to the introduction of the precise 
notations and assumptions needed in the sequel. In Section El a closed equation which governs 
populations only is derived from the Bloch equations up to an approximation of order 0(e 2< - 1 ~ M - > ). 
This equation is already a Boltzmann type equation, but the transition rates depend on time 
and £. The off-diagonal terms are proven to be negligible, of order 0(e l ~^). This transformation 
uses classical arguments for the Bloch equation in the weak coupling regime (see |Cas99| ICas02[ 
IHasOl] and also |KL57l IT7K581 lK7e83l l/wa66| for this point) . 

In Section HI a new Boltzmann type equation is derived removing the time-dependence of the 
transition rates up to an approximation of order 0{e 1 ~ 2 ^). These transition rates still include 
some terms which do not correspond to a resonance between the wave and the matter. These 
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contributions are removed in Section El with a new error of order 0(e^). The thus obtained 
equation is the rate equation we want to derive in this article. The main theorem in this Section 
(namely Theorem IHJ) yields the form of the transition rates in the rate equation. The techniques 
used are those of the averaging theory for ordinary differential equations (see [LM88, SV85 ). 
Diophantine estimates play naturally a key role in the analysis. Lemma [7| ensures that although 
the Diophantine condition is not stable with respect to small perturbations, violations of the 
condition only occur for large values of the indices, that are compensated by extra smoothness 
assumptions (see |CCC + 03] for a similar argument in another context). 

The sequel of the article addresses the analysis of the limit process e — > in this equation. 
Section El is devoted to the case study according to the value of fi/p leading to Tabled] in the 
finite dimensional case (finite number of levels) and to partial results in the infinite dimensional 
case (see Table EJ). The justification of the restriction to a finite number of levels is given in 
Section we show that there exists a number N such that considering only the N first levels 
implies an additional error of order o(l). The number N a priori depends on e, except when 
resonances do not occur beyond some level. 

Finally, we give in Section |S] the proofs of several lemmas, concerning continuity and non- 
positiveness of the relaxation operators first, then about existence and uniqueness of the as- 
sociated equilibrium state, and finally, we show the genericity of the Diophantine condition Q 
(see below). 



2 Functional setting 



The choice for functional spaces is first guided by some physical properties of the density matrix. 
Throughout this article the initial datum is taken such that 

p(0, n, m) = 0, if n ^ m, p(0,n,n)>0 and p(0, n, n) < oo. (4) 

n 

The energy levels are assumed to be bounded: 

Hn)) nm G r and (5(n)) nm G l°° 

(which is physically relevant, since these energies are bounded by the ionisation energy), as well 
as the relaxation coefficients 

( 7 (n, m)) N G inf 7 (n, m) := 7 > 0. 

We also suppose that the Pauli coefficients saisfy 

II^HWS^ := sup^|W(n,A;)| + sup^|W(n,A;)| < 00. 

k n . 

n k 

Matrix W describes the relaxation to a thermodynamic equilibrium. Entries W(n, m) and 
W(m,n) are therefore related by the standard microreversibility relation (see e.g. \BBR01\ 

W{n, m) = exp ( ' - )W(m, n), (5) 
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where T is a normalized temperature. This specific form is of great importance when describing 
the equilibrium states in Sections and 18.11 

We recall that the interaction matrix V is Hermitian. Thus we simply assume that 

llVlli-ji := sup^2\V{n,k)\ < oo. 

k 

n 

Classical ordinary differential equation arguments (see e.g. |Cas99j ) allow to state the 
existence and uniqueness of solutions to System (JTJ) for initial data in I 1 . Indeed, since we 
assumed that V(t, n, m) belongs to L°°(R+, n CO> and also that W(n, m) G H l^l^, 
we may use the estimate of Lemma The operators involved in Eq. (^J) are thus continuous 
on L°°(]R + ,Z 1 ). This implies in a straightforward way that these solutions exist and have the 
following regularity: p G C°(R + , I 1 ) and dtp G L°°(R + , I 1 ), for each e > 0. 

A weaker result is the existence and uniqueness in I 2 and we will have to restrict to this 
case from Section H] on. 

The summation and positiveness properties are preserved through the time evolution if the 
density matrix is solution to the Bloch equations More precisely we have the following 
lemma (points (ii), (Hi), (iv) are addressed in |Lin76|. lBBR01| ICasOlj ). 

Lemma 1. Let p(t = 0) satisfy conditions (jlj). Then, under the above assumptions, for all 
t E R, 

(i) there exists a unique solution p G C°(1R + ,/ 1 ) to Eq. (|T|). 
(ii) p(t) is Hermitian: p(t,n,m) = p(t,m,n)* , 
(Hi) the trace of p is conserved: 

pit* n ' n ) = n > n ) < 00 » ( 6 ) 



(iv) positiveness of populations is conserved: p(t,n,n) > 0. 

We stress the importance of items (Hi) and (iv), first established in |Lin76j . They give a 
nontrivial I 1 estimate for the diagonal part p&. This proves to be crucial in Proposition |21 (see 
also |Cas01j for a situation where the oscillations are much more difficult to handle). 

The electromagnetic wave is real and bounded: 4>(t) G R and <p L°°(R). Although a 
larger class of high frequency waves is addressed in [BFCDOB] , we have to restrict here to the 
quasiperiodic case. The amplitude is characterized via its frequency vector uj G R r and its 
Fourier coefficients <p a : 

(f)(t) = 4> a exp(ia ■ cot) where 0* = (j)- a and \4>a\ 2 < oo. (7) 
Here we denote a ■ uj :— a\toi + ■ • • + a r u r . 
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We need a certain number of assumptions on the wave and the interaction coefficients. 
Those already stressed will be valid (and not recalled) in the whole article. However, in some 
places, we will have to reinforce the decay assumptions and therefore we also label some specific 
hypotheses. 

As it is usual in the field of oscillations in ordinary differential equations and averaging 
techniques (see [Arn891 ISV851 ILM8 8 ) . we introduce a Diophantine condition on the frequency 
vector uj. 

Hypothesis 1 (Diophantine condition). There exists a (small) number rj > 0, and a con- 
stant C v > 0, such that 

Va = (oti, . . . , a r ) G Z r \ {0}, V(n, k) G N 2 such that a ■ uj + uj(n, k) ^ 0, 

C 

\a ■ uj + uj(n, k)\ > - ; — - — — — ; — - — — - — ; — — , (8. a) 

and similarly 

VaGZ r \{0}, \a-u\>- °l - . (8.b) 

i (1 + |a|) r - 1 +'' 

Remark. Given once and for all a fixed 1] > 0, we can classically claim (see I Am Sty ) that 
there exists a constant C v > 0, depending on uj(n, m) and on rj, such that for almost all value of 
the frequency vector uj = (ui, . . . , u r ) Hypothesis^ is satisfied. This is proved in Appendix \8.Z& 
This condition is therefore not much restrictive. 

The same kind of estimate is also needed in the case when a = 0, which means that the 
energies do not converge too fast towards the ionisation energy. 

Hypothesis 2 (Convergence towards the ionisation energy). There exists a constant 
C v > 0, such that 

V(n,A;)GN 2 , \u(n,k)\>- — — — , 

where rj is the number occurring in the Diophantine assumption (Hypothesis^. 

From Section 0] on, we impose the following hypothesis in order to average transition rates 
in time. 

Hypothesis 3 (Smoothness assumption). The Fourier coefficients 4> a satisfy 

^(l + |a|r 1+ l0 a | 2 <oo, 

a 

where rj is the number occurring in the Diophantine assumption (Hypothesis^. 

From Section El on, we need the following two hypotheses to justify the sorting out of 
resonant and non resonant contributions. 
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Hypothesis 4 (Reinforced smoothness assumption). There exists N v > 2p/p such that 
the Fourier coefficients <p a satisfy 

^(l + |a|)^ 1+ ^|0 Q | 2 <oo, 

a 

where rj is the number occurring in the Diophantine assumption (Hypothesis^). 
Hypothesis 5 ("Far from continuum" assumption). The interaction coefficients satisfy 

supj^ ((l + n) 1+ "(l + m) 1 - H ') jV " \V(n,m)\ 2 < oo, 

n 

m 

where rj is the number occurring in the Diophantine assumption (Hypothesis^) and N v is given 
by Hypothesis^ 

This last hypothesis means that only low levels (i.e. levels which are far from the continuous 
spectrum or ionisation threshold) really interact with the wave with a significant contribution. 
To restrict to a finite number of levels we also use an assumption on the interaction of low and 
high levels via the relaxation operator. 

Hypothesis 6 (Weak interaction of low and high energy levels). The longitudinal 
relaxation coefficients satisfy 

sup^(l + n) K [l +m) K \W(n,m)\ < oo, 

n 

m 

for some K > 0. 

We introduce a last notation useful for describing approximations of p. 

Notation 2. For s > 0, and q > 1, the symbol Ol™ qr,z«) {s s ) means that for all T > 0, there 
exists a constant C > 0, that does not depend on e, such that the corresponding term is bounded: 

\\Ol°°({o,t],ii)(£ s )\\l°°(Io,t],i<i) < Ce s . 

In the whole article C denotes constants which do not depend on the (small) parameter e. 
It however possibly depends on all the coefficients of the problem and on the initial data, but 
we will never make this dependence explicit. 

3 A first closed equation for populations 
3.1 Populations and coherences 

In the same way as we defined the populations p^(t,n), we define the coherences as the off- 
diagonal part of the density matrix 

Pod{t, n i m ) — p(t, n i Tn)l[n 7^ m]. 
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Set 

Q e (n, m) := —iu(n, m) — ie p 5(n, m) — e^ r y(n ) m), 
where we stress that Q £ (n, n) = 0. With these notations, Eq. JIJ reads for the populations: 



V ( -z,n,k J Pod(t,k,n) - V I -z,k,n ) p od (t,n,k) 



k 

+ ^yV{k,n)p A {t,k) -W{n,k)p A {t,n)\ 



(9.a) 



and for the coherences: 



dtPod{t,n,m) = —Q £ (n,m)p od (t,n,m) 
e z 

+ -V ( -2,n,m ) [p d (t,m) - pd(t,n)] 



(9.b) 



e \e 
i 

e 



V ( — ,n, fc J p od (t,k,m) - V (\,k,m ) p od (t,n,k) 



As a consequence of the Hermitian properties of Eq. (0) recalled in Lemma ^ Eq. ()9.a|) can 
also be cast as 



d t pd(t,n) 



-- Im 



V [—,n,k) p od (t,k,n 



(10) 



Y^[W{k, n) P d(t, k) - W{n, k)p d (t, n)). 



3.2 An equation for populations only 

In this section, we transform the coupled system ()9.aj) - (J9.bj) into one equation governing the 
populations p d (t,n) only. More precisely, we show the following proposition. 



Proposition 2. Define the time dependent transition rate 

-t/e 



V £ (^k,nj :=2\V(n,k)\ 2 Re dsexp({l e 



(k,n)s)(f)[ — )(/)[— -s 



Then, for all T > 0, the vector p d satisfies 

d tPd (t) = (* e + W^j p d {t) + L ~ i[W i) {e l ~ 2 ^) ■ (11) 

Remark. Eq. (|llj) is a linear Boltzmann type equation with a time dependent transition rate. 
This is our first description of the behavior of the populations at leading order in e. 

We could likewise obtain approximations of p d at each order (in e), thus providing a hier- 
archy of Boltzmann type equations for the successive orders. 
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Remark. Using Lemma [771 and since KeQ £ (n,m) < — £ M 7, uniformly in n and m, n 7^ m, 
operator \l/ £ is a priori of order e^ 1 - 1 on I 2 or I 1 . More precisely we have the following estimates. 

Peit/e^ml, < Ce-^||y(n,A:)||?„,i nJ « )I x||0||iac CR) ||pd||ti, (12.a) 
\\%(t/e% Pd \\ l2 < Ce-^ll^^llf^i^ll^llloo^llpdllp. (12.b) 

Proof. Proposition's proved in three steps that follow BFCD03]. First, coherences p Q d(t, n, k) 
are computed at leading order in e. These leading order terms are expressed in terms of 
populations p d (t,n) only as stated in Lemma 01 below. Next, we plug this result into Eq. (jIDJ) 
governing populations. The closed equation for populations (jUj) which is thus obtained is a 
linear Boltzmann equation with a time-delay term. To show this delay is indeed small and 
can be removed to obtain the delay-free equation (fTT|) we then use calculations inspired by 
Cas99[IHa^02llCas01| . 



First step: computation of coherences 

Since the initial data for coherences is p d(t = 0, n, m) = 0, the integral form for Eq. (|9.b|) reads 



ft/e 2 

p od (t,n,m) — ie dsexp(Q, e (n : m)s) 
Jo 

+ (A £ p od )(t,n,m), 

where 

t/e 2 



V - s, n, rn^j [p d (t - e 2 s, m) - p d (t - e 2 s, n)) 

(13) 



(A £ p od )(t, n, m) := ie I ds exp (Q £ (n,m)s) x 



x 



E 



V ( — - s, n, k ) p od (t - e 2 s, k,m) - V ( -= - s,k,m) p od (t - e 2 s, n, k) 



We are only interested in the first term in the expansion of p od . Indeed, we only use leading 
order terms in e. From the integral equation (JHH), we would however be able to express p od 
as a complete expansion in powers of e, in terms of p d . The following lemma states that the 
remainder of the expansion is indeed small. 



Lemma 3. Set 



Pol {t,n,m) := ie J dsexp (Q £ (n,m)s) V ^— — s,n,mj [p d (t — e 2 s,m) — p d (t — e 2 s,n)}. 

Then, for any given time T > 0, there exists a constant C , that does not depend on e, such 
that we have the estimates 



t/e 2 



(0) 

Pod - P od 



Pod 



Z,°°([0,T],U) 



L^dC^T],; 1 ) 



< Ce 2 ^\ 



< Ce 1 -^. 
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Proof. (Lemma El)- 

Let T > be given. Because Ref2 £ (n, m) < —£^7 < 0, uniformly in n and m, n ^ m, we have 
the estimate 



ds |exp (Q e (n, m)s) 



o 

1-fMl 



L-(R+,C4n/S?ii)IIPod||L-([0,T],Zi) 



||^4ePod||L°°([0,T],/ 1 ) < 2£ 

< Ce L At ||pod|U°=([o,T],« 1 ) 
According to the definition of p ° d \t,n,m), Eq. (fT^|) 

Pod (t, n, m) = p °/ (t, n, m) + ( A £ p od ) (t, n, m) , 
and we can estimate the difference 

1 1 Pod — P d IU°°([0,T],U) < At ||Pod||L°°([0,T],Z 1 ) 

< Ce 1 ~^ 1 1 p od - p o ^ ||i«>([o,T] ) i 1 ) + Ce^llpidlU^IO.Ty 1 )- 



Hence, if e is small enough, we have 

II CO) 1 1 ^ n 1— mm CO) || 

1 1 Pod - Pod ||i oo ([0,T] ) Z 1 ) S llPod IU oo ([0,T] ) Z 1 )- 

According to the definition of p o d \t,n,m), and using the same type of estimates, 



l^id IU°°([0,Tyi) < Ce l M ||V|| iO o( K;i oo;i n ;o Ci i)||pd||Loo([ 0)T ] i/ i). 

Now the crucial estimate stems from the trace conservation property © which also reads 
||Pd||.L° c ([o,Ty 1 ) = ||Pd(£ = 0)||ji. This / 1 -estimate on the diagonal terms of the density matrix is 
therefore sufficient to control all the off-diagonal terms in turn, and Lemma El follows. □ 

Second step: the time- delayed differential equation for the populations 
Lemma H3 together with Eq. (JTUJl governing pa implies that 



d tPd (t,n) = ^yV{k,n)p A {t,k)-W{n,k)p A {t,n)} 

k 

f t/£2 

+ 2_. / ds[p d (t — e 2 s, k) — p d (t — e 2 s,n)] x 

r. JO 



(14) 



x2 Re jexp(ff(£;, n)s)V n, k^j V — s,k,n 



+0 



L°°([0,T],Zl) 
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Third step: convergence to a delay-free equation 

From now on we will only deal with Boltzmann type equations and use extensively the shorter 
expressions defined in Notation ^ Hence Eq. can be cast as 

d t p d (t,n) = (^ £ {t/e\p d ) (t,n) + (W s p d )(t,n) + L ^ {[oW (e 1 " 2 ") . 

Moreover, if T > is given, the delayed terms p d (t — e 2 s) in Eq. (fTljl read 

p d (t - e 2 s,n) = p d {t,n) + O (er 2 s||5tp d (-, rz) |U°°([o,^])) • 

Thus, Eq. (PJ) yields 

d t Pd(t,n) = (^ £ {t/e\p d ) (t,n) + (W iPd )(t, n) + 0^ ([(W) (e 1 " 2 ^) + r £ (t, n), 

where the remainder r e can be estimated by 



kelU^O.T],/!) < C||ftPd||L°°([0,T],Z 1 ) 



< Ce 2 -^e 

< Ce 2 -^ 



tPd||L°°([0,T],/ 1 ) 
|Pd||L°°([0,T],Z 1 ) 



ds s exp(f2 £ (n, m)s) 



thanks to Eqs JIEJ) and |T2~aj) 

thanks to Eq. (0), 



(15) 



for some constant C, that does not depend on e, if e is small enough. Including this new 
estimate in Eq. (fT5]l. we have 



d t Pd(t, n) = ((Ug tt p d ) {t, n) + (W iPd )(t, n) + O l ^ {[W) (e 1 " 2 " + e 2 - 3 ") 
We recall that p < 1/2 therefore 1 — < 2 — 3/x and Proposition El follows. 



□ 



4 Time-averaging of transition rates 

Proposition |2] reduces the problem to the asymptotic analysis of a leading order equation, 
namely Eq. (fTTj) . As e goes to zero, the rapid oscillations of the time-dependent coefficient 
^ £ (t/e 2 ) are naturally smeared out, so that only the time average of ^ £ plays a significant role. 
The corresponding rigorous mathematical statement is proven here, as in [BFCD03J, upon using 
averaging techniques (see e.g. [SV85J). This leads to substitute Eq. (fTTj) by a new Boltzmann 
equation, with time-independent rates. 

The result of these averaging procedures strongly depends on the specific form of the wave. 
In the present article, we need explicit forms for the averaged transition rates and therefore 
restrict to the case of a quasiperiodic wave, as described in Section El The small divisor 
problems that stem from this analysis are handled assuming strong decay properties on the 
Fourier coefficients, namely Hypothesis 01 

The averaged transitions rates we obtain in this section still depend on e. Besides we do 
not discriminate dominant and negligible (with respect to e) contributions. We postpone this 
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sorting out to Section |3] where negligible contributions will be pointed out as non-resonant 
contributions. 



Using the explicit value of the wave 4>{t) given by Eq. (JJJ), we compute from Eq. (JTTjl the 
time-dependent transition rate \l/ e : 



tyJ — ,k,n\=2\V(n,k)\ 2 Re ^ Mpexp (i(a + p) ■ u—j x 



1 - exp ([-e"7(fc, n) - ^g, n) + g ■ + eH{k, n))]t/e 2 ) 
X e^7(fc,n) + i(o;(A;,n)+/?-^ + £: p 5(fc,n)) ' ^ ' 

The goal of this section is to prove the following result. 
Proposition 4. Define the averaged transition rate 

i r 

(^ e )(k,n) := lim - / ds\l/ e (s). 

S-++OC, s Jo 

Its explicit value is 

(%)(k,n) = 2\V(n,k)\ 2 Y - - , ^ ( , fc,W l r/ , ,,J ^| 2 . (17) 

X /3 ^.^ 7 (A;,n) 2 + ^(n,A;)+/3-^ + ^(fc,n)| 21 ^ 1 

Define as the vector solution to 

dtti ] = «*.) + W)„pf } (18) 

wift initial data p^\o,n) = p^(0,n). Under Hypotheses Q] and\Q for all T > 0, there exists 
C > stzc/i t/iat 

II Pd — /^IU°°([0,T],J 2 ) — C^ 1-2 ^- 

Remark. In the transition rate (\& e )(fc, n) there are still contributions of non-resonant waves, 
i.e. thoses such that ui(n, k) + /3 • u ^ 0. JTiq/ are removed in the next section. However, we 
gained in defining a time-indepedent transition rate, in that we filtered oscillations. 

The simple but crucial remark that leads from Proposition El to Proposition 0] is the fact 
that (* e ) e V£l\ n Ifli and its entries are non-negative. We can therefore apply Lemma ITT1 
directly. 

Lemma 5. For any e > 0, operator {$ £ )$ is a bounded non-positive operator on the Hilbert 
space I 2 . In particular, the exponential exp(t(^ e ) j) is well defined as an operator on I 2 fort > 0, 
and its norm is less than 1, for all t > 0. 

(Proof: see Appendix 18. 11 Lemma fTT|) 

Remark. This property proves to be crucial in the sequel, and therefore the asymptotic result 
in Proposition ^] and Theorem may not hold when the transition rate has no sign. This is 
why we restrict the analysis to the case of a quasiperiodic wave when relaxations tend to zero 
with e. The asymptotic result may also be false if applied to the time- dependent transition rate 
^ £ (t/e 2 ), which clearly changes signs. 

Proof. We follow [BFCD03 to prove Proposition 0J 
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First step: splitting of \I/ e 

Proposition 0] amounts to estimating the difference 

A(t) :=p d (t)-p^(t). 

To this aim, we first give shorter forms for the equations governing p d and respectively. 
Namely Eq. (fTTj) can be cast as 

dtp d (t) = U e (j^jA (*) + (^Pi)(t) + {e 1 - 2 ^) . 
where the transition rate ^ £ {t/e 2 ) is given by Eq. (|16|h and Eq. ()18|) also reads 

Hence the difference A(i) satisfies the equation 

$A(f) = «*.),A)(t) + ^ sc (Ji)^) (*) + (Wj|A)(t) + 0(e 1 -^), (19) 

where 

contains the oscillatory contribution to the transition rate, which we want to prove to be 
negligible. Gathering the terms for which a + (3 = 0, this contribution is equal to 



^(^k,nj = 2\V(n,k)\ 2 Re( -J2 



, r £^ 7 (A;, n) + i(u(k, n) + (3 ■ u + eP5(k, n)) 

x exp([-e M 7(fc,n) - i(uj(k,n) + (3 ■ lu + e p 5(k,n))}-!-) 
a 0/3exp(i(a + /j) -^t/e 2 ) ^ 

x [1 - exp([-e M 7(fc,ra) - i(u(k,n) + (3 ■ to + e p 5(k,n))]^)] 

This expression carries "time-oscillations" (at frequency e~ 2+fl at least), which kill the possibly 
diverging factors (due to the denominators), and make them of size e 2 ~ 2 ^. 

Second step: preliminary bounds 

Since V G n and ^ | <^>/3 1 2 < oo, we first find that 

sup£K¥ e )(fc,n)| <Ce~^ sup ^ |*° sc (t/e 2 ,M)l < Ce'", 

k k 

sup V|(f e )(fc,n)| <C£-^ supV|^ sc (t/£ 2 ,fc,n)| < Ce - ", 

n n 
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for some C > that does not depend on t and e. Lemma ITT1 yields the operator estimates 



Besides, we have the upper bound 



£ H 



P 



< Ce~nu\ 



(21) 



sup 

0<t<T 



t/e 2 



ds^° e sc (s,k,n) 



< Ce 



-2fi 



(22) 



l n ''k li k n 



Now, according to Eq. (|20|). \I/° sc (t, fe, n) is a sum of two different terms. We use the decay 
assumptions V G fl /j^Z* and a G / 2 to estimate the contribution of the first term to the 
integral by 



C 



\e^(k, n) + i(oj(k, n) + (5 ■ u + e p 5(k, n)\ 2 
The second contribution is estimated by 

\V{n,k)\ 2 \U\h\ 



< Ce-^ L . 



/oo/lr-i/ 00 /! 
ri L k lL k ri 



C 



\e»y(k,n) + i(u(k,n) + (3 ■ lu + eP6(k,n))\ \(a + 0) ■ u\ 



/oo/lo/oo/l 

n fc 1 fe n 



-c 



\V{n,k)\ 2 \<i>o\\M 



\e^(k,n) +i((a + f3) ■ u + eP5(k, n))\ : 



n l k' 1 k n 



a, f3 

< Ce~ 2 » 



\a 



a,P 



thanks to the Diophantine estimate (Hypothesis [IJ, together with Hypothesis 01 This yields 
inequality 



Third step: integral form of the equations 

Since A(0) = 0, the integral form for Eq. (fTT?)) governing A(t) reads 



where 



A(t) := A (1) (t) + A (2) (t), 

Jo V£ y tt 

A^(t) = r^exp([t- s ](vI/ £ ))((H/ fl A)( S ) + 0(£ 1 - 2 ^)). 
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Fourth step: estimating A^'(t) and A^ 2 '(t) 

Here Lemma |5] proves crucial, in that we use the exponential of the bounded operator (^ e ). 
This together with Lemma ITT1 applied to operator W yields an estimate for A^: 



\\A®(t)\\ P <C[e 1 -*»+ / ds\\A(s)\\ P . 



(23) 



On the other hand, to take advantage of the time oscillations of the operator ty° sc (t/e 2 ), we 
carry out a natural integration by parts in the expression for A^: 



A«(t) 



t/e 2 



du*r(u) paw 



+ e 2 J dsexv([t-s]{* £ )){* e ) I J »° SC ( M ) j ,,,,(*) 

- e 2 dsexp([*-s](tt B )) (j^ »° sc («)j x 

x (<tf e ) + M/° sc (A) + W + Oie 1 -^)) PA (s), 

where we have used Eq. (JTTJ) to express d t pd(s). Estimates (J2T|) on the operators (\l/ e ) and 
\l/° sc (t) together with Lemma ITT1 and Lemma El (non-positiveness of ($? £ )) lead to 



\\A^\\ Laa{mP) <Ce 2 ^ sup 



0<t<T 



t/e 2 



||Pd|U°°([0,T],/ 2 )- 



C(l 2 ) 



Besides we have ||pd||L°°([o,T],z 2 ) < ||Pd||L°°([o,T],« 1 ) < C, and it follows that, for e small enough, 
we have 



|A (1) || L oo ([0>T]j/2) < Ce 2 -» sup 



0<t<T 

< Ce 2 -» sup 

0<t<T 



t/e 2 



dsV° sc (s) 



C(l 2 ) 



t/e 2 



ds^° £ sc (s } k,n) 



t 7l ''k' ' L k n 



< Ce 2 -'^. 

This, together with estimate (J2BJ), and Gronwall lemma, yields 

\\A{t)\\ L ~ {[W) <Ce l - 2 ». 

and Proposition 0] is proved. 



□ 
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5 Keeping only resonant contributions in transition 
rates: the main theorem 



In this section we prove that the non-resonant contributions, which correspond to the triples 
(n,k,/3) such that u(n,k) + [3 ■ u> ^ in the transition rate (fTTjl . are negligible in the limit 
e — > 0. We therefore replace the transition rate (\& e ) by a purely resonant transition rate 
(\l/ e ) dom , which however still depends on e. To get rid of this last dependence, we will have to 
specify the value of p/p, which is the goal of the next section. 

Still due to small denominator problems, we need to reinforce the decay assumptions on the 
coefficients and assume Hypothesis 0] and |S] hold. 

Theorem 6. Define the transition rate 

,v ' y v ' ' P&Z r ; w(fc,n)+/3-w=0 

(2) 

Let a/so fre solution to 

dtp? = ((**)*** + W) t p® (25) 

itra£/i initial data p d 2 \o,n) = p^(0,n). We assume that p < 1/2. Then, under HypothesesUl\R 
^] and /or a// T > 0, i/iere exists C > swc/j t/iat 

llPd - IU<-([o,n,P) < C(e» + e 1 -^)- 

Remark When p = 0, Theorem^ does not give a good approximation of pd, and we do not have 
a better description than the one from Proposition^ where all (resonant and non-resonant) 
frequencies have a contribution to the transition rates. 

The proof of Theorem El follows closely that of Proposition However, we notice that the 
resonant values that play a role here satisfy u>(n, k)+/3-uj = and not u{n, k)+j3-ui+e p 5(n, k) = 
0. Therefore we need to understand the effect of a perturbation on Diophantine estimates. 
Lemma [7| below answers this problem. 



5.1 Perturbed Diophantine estimates 

Lemma 7. If u and u(n,k) satisfy the Diophantine condition (18. a|) and Hypothesis \B with 
constants r\ and C v , then the following assertion holds. 
If (n,^)eNxNx Z d satisfies 

\f3-LU + uj(n,k) +e p 5(n,k)\ < 1 " 



2 (1 + \(3\y- 1+ ^(l + n) 1+ v(l + fc)i+r 
then 

(1 + \{3\y- 1+v (l + n) 1+ \l + k)^ > 
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Remark. The Diophantine condition (|8.a|) is not stable with respect to small perturbations: 
coefficients u(n,k) + e p 5(n, k) , that can be arbitrarily close to u(n,k), are capable of violating 
the Diophantine condition (J8.aj) . Indeed, for e small, the condition 

\f3-u + u(n,k)+e p 5(n,k)\ > 1 " 



2 (1 + \(3\) r - 1+ v(l + n y+v(l + k) l+r > 



does therefore not necessarily hold, only assuming that the left hand-side is non-zero. Neverthe- 
less, Lemma^ claims that this condition may only be violated for values of the triple (n,k,(3) 
which are very large when e — > 0. See \CCC*~ 05lf for a similar argument. 

Proof. Set K = C v /2 and take (n, k,/3) such that 
\P-u + u(n,k) + e p 5(n,k)\ < 



(1 + \{3\y- 1+r ](l + n) 1+r >(l + ky+v 
Then 

\(3- u + u{n,k)\ -e p \5(n,k)\ < h 



(1 + \f3\y- 1+1 i(l + n y+v(l + k) 1+ n : 
and according to condition (|8.aj) (or Hypothesis when (3 = 0) 

2h e p \5(n,k)\< R 



(l + lPiy-i+n^ + ny+v^ + ky+n 1 v ' n - (i + \p\y- i -*>(i+n) i+r i(i + k) i -*>' 

Hence 

h <e p \5(n,k)\, 



(i + \(3\y- i+ v(i + n y+v(i + ky +r > 

which ends the proof of Lemma □ 

5.2 Proof of the main theorem 

First step: Integral formulation 

To deduce Theorem |H1 from Proposition 0] only amounts to estimate the difference 

A(t)=pi 2) (t)-pf ) (t), 

where we use once again the notation A. Now pl\t) and p£\t) are respectively solution to 

dt P f{t) = ((^)f om pi 2) ) (t) + (W,pf) (t), 

and 



Hence 



d t A(t) = ((%hA)(t) + ((^V?) (t) + (Wj,A)(f). (26) 
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where 

(* £ r g (M) = (Ve) dom (k,n) - (* e )(k,n) 

contains the contributions to the transition rate, that we want to prove to be negligible. Since 
A(0) = the integral form for 1)2 6 j) reads 

A(f) =jfdsexp([t -*]<*.)) ((^)r g pf) (s)+iy S eM[t-s](%})(W i A)(s). 

Second step: Estimating (\E' £ ) ncg 

In view of Eqs (JTTj) and (|2^jl , we have 



2\V(n, k)\ 2 e^j(k,n)\cP P \' 



" m ^e^7(M) 2 + MM)+/^ + e p *(M)| 2 



We split this expression into two contributions according to the fact that 



\P-w + u{n,k) + e p 5(n,k)\ > -- 



2 (1 + ^^-^(l + + jfe) 1 "^' 

or not. Using Lemma [7| for the second contribution, we obtain 



/oo/lri/ooH 
' n b U "fc 'n 



< 



sup{ y: mn ' k)?£ c {K n) 10/3,2 a + iflr 1 ^ + + fc ) 



fc,/?; w(n,fc)+p'-u;^0 

+ £ 1 

fc„3; w(n,Jfc)+/J-u^0 



(1 + \(3\) r - 1+ri (l + n) l+ri {l + k) 1+r > > Ce~ p 



l+ri 



2\V(n,k)\ 2 £^\<j> p \ 
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The first sum is estimated using Hypotheses EJ 121 and E| The second term is first multiplied 
and divided by the quantity [(1 + \(3\) r - 1+ ' n (l + n) 1+r *(l + k) 1+r >) Nr >. Therefore we get 



n "k' ' k n 



n,k,P 



r— 



(l + n^l + jfe) 1 -^)^ \V(n,k)\ 



2 1 / |2 



Since we assumed that > 2/i/p, we finally have 

\\(%r%^ llniri i <Ce». 



Third step: Conclusion 

Using Lemma El (non-positiveness of (\& e )) leads to 



IA 



< C\\(^ £ ) ae ^c(P)\\pf ) \\L^([onP) + C f ds\\A(s)\\ P 

Jo 

< C\\(^ £ ) neg (n, k)\\i^ii n i^ii || A || loo([ ,t],« 2 ) + llPd" IU°°([o,T],i2) 



C / ds||A(s) 
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A by-product of Proposition 0] is that the quantity ||p d ||l^([o,tu 2 ) can be estimated by a 
constant. Therefore 

||A(*)||, a < C|K* e )^(n,A:)||^i nl -,i(l + ||A|| L =o ([ o, nia) ) + (7 [ ds\\A(s)y 

Jo 

< Ce^l + \\A\\ LOO{W) ) + C [ ds\\A(s)\\ l2 . 

Jo 

Thanks to the Gronwall lemma 

||A||iac([ 0j T],/2) < 

which proves Theorem |3 □ 

6 Time-layers and equilibrium states in the e —> limit 

In the previous sections we have derived the rate equation (|25jh with the transition rates 

W™ d (n,m) = (q £ ) dom {n,m) + W{n,m), 

which may be considered as modified rates, via the interaction with the wave. We now turn to 
the study of the dynamics of the solution to Eq. (J2o]) . 

For a fixed e, this dynamics is described in Appendix 18.11 However, in most cases, the 
coefficients prove to be singular in e. Therefore the time evolution of the solution obeys two 
different regimes as e goes to zero: first a time-layer, and then relaxation to an "equilibrium" 
(if it exists, see Appendix 18.1)1 . The duration of the time-layer is always less than 0(1), thus 
the Bloch equation and the rate equation (|25[) behave in the same way on this time-range 
(as Theorem El asserts). On the contrary, the relaxation towards an equilibrium state, as time 
goes to infinity, is a priori specific for the rate equation, since Theorem El only applies on the 
fixed time-interval [0, T]. 

In this section, we split the rate operator (W™ A )§ into three contributions. 

1. Some levels are decoupled from all other levels. The corresponding columns (and subse- 
quently lines) in the matrix (W / e mod )(j are identically zero. These levels have a constant 
population for all time (exactly for the rate equation, and at leading order for the Bloch 
equation). They are excluded from the sequel of the argument. We call E the subspace 
spanned by these levels in I 2 (this subspace does not depend on e). The projection onto the 
levels which have a non trivial time evolution is denoted by Ilo. Therefore So = Kerllo. 

2. We split (\l/ e ) dom into two contributions: (\l/ e ) smg collects the singular coefficients, which 
go to infinity as e — > 0. Then, (\[/ £ } nonsm g gathers the non-singular coefficients, which 
are 0(1). The corresponding case study according to the value of fi/p is performed in 
Section 16.11 

3. Finally we define the / 2 -orthogonal projection II onto the space Z 2 -orthogonal to S in 
Ker(\l/ £ ) smg . It does not depend on e either. The quantity (1 — II) pa will be proved 
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to vanish in time 0(e a ), where a depends on the ratio p/p (see Table [Q). After this 
time, i.e. past the initial time-layer, a polarized solution p^ Dl persists. It is given by the 
£-independent system 

u P f = pf, d tP f = uwf'upf, 

where W?° l = W + ^ onsin s. 

This section is organized as follows: in Section 16.11 we set some notations and explain why 
the transition rates are in general splitted into three types of terms. Their respective sizes are 
given by powers of e. Under the assumption that the number N of quantum levels is finite, a 
precise description of the time layer is given according to the value of p/p in Section RHfl The 
key lemma for this analysis is given in Section 16.21 Section 16.41 is devoted to the particular 
case p = 0. In Section l63| we discuss the few cases (in terms of values of p/p) when we may 
conclude for an infinite number of levels. A rigorous framework for the justification of the 
restriction to a finite number of levels is given in Section 



6.1 Setting for the time-layer 

We first split (\E f e ) dom (n, m) into two contributions according to whether 5(n,m) = or not. 
We restrict the discussion to the case when p > since the form of the transition rates in 
Section does not apply when p = 0. This case is treated separately in Section 16.41 To 
simplify notations, we set 

C(n,m) = 2\V(n,m)\ 2 ]T \<f> \ 2 , 

f3€Z r ; u>(n,m)+f3-u>=0 

and (\l/ £ ) dom (n, m) = A e (n,m) + B e (n,m), where 

A E {n,m) := C(n,m) - - ^jfa^ ^l(S(n,m) = 0), (27.a) 

B B (n,m) := C{n,m) _ . ^^'^l, ^l{5{n,m) ^ 0). (27.b) 

Of course Eq. (|27.a|) also reads 

AJn, m) = e~^C(n, m)—^- — rl(8(n, m) = 0) =: £ _/ M(n, m). 

j(n,m) 

In a similar way we rewrite Eq. (j27.b|) as 

B E (n,m) =: e~ lJ B £ (n,m), 

where v is chosen such that B £ — > B° in I 00 ! 1 , as e — > 0, and therefore also as an operator on 
I 2 . To this aim we consider two cases, namely p < p and p > p. We set 

v = p and B £ (n,m) := C{n,m)— — — ' -l(5{n,m) ^ 0) if p < p, 

7(n, my + e 2[ P Wd)(n, m) z 

v = 2p - p and 5 e (n, m) := C(n, m) ' m) / 0) if p > p. 
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With the above notations, we have cast Eq. ()25j) governing y = p^' as 

d t y = (e-»A + e-»B e + W)0. (28) 

Using the block decomposition described in Appendix 18.11 we notice that Ker(B e ) = Ker(i?°) 
is constant with respect to e > 0. This is a priori only true if all the minimal stable eigenspaces 
are finite-dimensional. 

Now, sorting between singular and non-singular contribution yields the different cases in 
Tabled In the framework of BFCD03J B £ is identically zero. If p < p then B e contributes to 



singular terms at the same order as A. If p < p < 2p, there are two orders of magnitude 
and e~ v in the singular term. Finally if 2p < p, the effect of B e is of order 0(1) and 
must be associated to that of W . In other words, if p < 2p, the non-zero entries in A and B e 
contribute to (\l/ e ) smg , and (v|/ e } nonsm g = 0. If p > 2p, only the non-zero entries in A contribute 
to (\l/ e ) smg : (v]/ £ ) nonsm s = B £ , which has a non-vanishing value \j/jij onhmg only if p = 2p. 



6.2 A finite dimensional lemma 

Lemma 8. Let < v < p. Let A and B £ G .MjvOR) be symmetric non-positive matrices 
such that B £ — » B° , and assume that Ker(B £ ) is constant for e > 0. Let U be the orthogonal 
projection onto Ker A fl Keri? . Then there exists a constant c > 0, namely 

c= min —((A + B°)x,x), 
\\x\\<i,nx=o 

such that any non-zero eigenvalue X s of (1 — T\){e~^A + e~ u B £ ){l — II) satisfies 

\ £ < -ce- v . 

Proof. Since p > v, we write 

(1 - n)(e"M + e- y B){l -U)= ^(1 -U)(B + e^A^l - II). 

Let k £ be a non-zero eigenvalue of (1 — U)(B £ + e u ^^ A){1 — II). Since A and B £ are non-positive 
matrices, we have k £ < 0. We therefore want to prove that for the constant c defined in the 
lemma, k £ < —c. Let x £ G Range(l — II) be an eigenvector associated with the eigenvalue k £ 
and such that ||a; £ || = 1. We have (1 — II) (B £ + e u ~^A)(1 — Il)a; £ = k £ x £ . Taking the scalar 
product with x £ , we obtain 

k £ \\x £ \\ 2 = {{B £ + e u ~^A)x £ ,x £ ) = (B £ x £ ,x £ )+e u -»(Ax £ ,x £ ). 

Since e G [0, 1], v — p < and (Ax £ ,x £ ) < 0, we have e u ~^(Ax £ , x £ ) < (Ax £ ,x £ ). A subsequence 
of x £ converges to x° and therefore 

sup k £ \\x £ \\ 2 < ((A + B^x^x ) < -c\\x°\\ 2 . 

£>0 

Finally, if c were zero, then ((A + B°)x,x) = that is (Ax, x) = and (B°x, x) = and 
therefore x G Ker A fl Ker B°. Thus (1 — II)a; = and x = which is impossible. □ 
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Remark. 

(i) The above proof is restricted to the finite dimensional case for two reasons. First, in the 
infinite dimensional case the convergence of x £ to x° would only be weak. Second, the 
maximum of ((A + B°)x,x) on vectors x such that \\x\\ = 1 could be zero. 

(ii) To avoid such limitations we could think about replacing B £ by the leading order terms 
(those which are not vanishing as e — > 0) in the series expansion of B £ (Lemma\^ can 
clearly be extended to a finite number of matrices) extending the computations of Sec- 
tion \6.FA to sums of several powers of e. This is not possible since such a procedure would 
not ensure the non-positiveness of the resulting operators, which is crucial in the proof. 

Lemma |H1 allows to split the solution y to Eq. (j2EJ) into two parts: (1 — H)y vanishes 
exponentially and a solution z associated with the initial data Hy(0) survives. 

Corollary 9. Let < v < [i. Let A,B e e Aipf(M.) be two symmetric non-positive matrices 
such that B 6 —>■ B° , and assume that Ker(B 6 ) is constant for e > 0. Let W G .Mjv(K), and let 
IT be the orthogonal projection onto KerAjj H Ker B^. If y is solution to 

d t y = (e^A + e~ v B £ + W\y 

and if z is solution to 

d t z = riWtflL?, z(0) = ny(0), (29) 

then 

||(l-n)y|| < C[e v + exp(-cte- u )), (30.a) 

My-z)\\ < C(e v + eM-cte- v )), (30.b) 

where C depends on W , 2/(0) ; T and c (the constant of the previous lemma). In particular, 
lim £ _ (l - U)y = 0. 

Proof. The proof is standard. We reproduce it for the sake of completeness. We first prove 
estimate (|3U.a|) . Since II is the orthogonal projection onto Ker A (1 KerB £ , we have ITAjj = 
nS| = A fl n = B £ U = 0, and therefore 

d t (l-U)y=(l~U)d t y = (l-U)(e-^A + 6- v B e ) iy + (l-U)W0 

= (e-^A + e^B% y +{l-Yi)W iy , 
= ( £ -^A + e^B £ Ul-U)y + (l-U)W iy . 

The solution to this equation reads 

(1-U)y(t) = exp(t(e-M + ^ J B £ ) fl )(l-n) 2 /(0) 

+ f dsexp((t-s)(e^A + e~ u B £ ) t )(l~U)W iy (s). 
Jo 
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Lemma |H1 yields a bound for the spectrum of (e fJ, A + e u B e )$, namely Sp(e ^ A + e u B e )§ < 
—ce~ v and we also have ||1 — n|| < 1, therefore 

||(1 - U)y\\ <eM-cte~ v ) \\y(0) \\+ C f dsex V (-(t - sjce""). 

Jo 

Finally ||(1 - U)y\\ < C{e v + exp(-cte~ u )) which is estimate (IHlTal) . 

We now prove estimate ()30.b|) . We write y — z = U(y — z) + (1 — H)y — (1 — H)z. By 
definition Hz = z and we already have estimated (1 — IT)?/. Thus there remains to estimate the 
quantity H(y — z), which is solution to 

d t Tl(y - z) = d t Tly - d t Uz = Ud t y - d t z 

= U(e-^A + e~ u B\y + UW iV - LW fl ILz 

= UW t y-UW t Ilz 

= IlW t Uy + UW t (l -U)y-UWflz 

= ni^n u( y - z) + nw$(i - U)y. 

The solution to this equation may be estimated by 

||I%-*)|| < [ ds\\ex V ((t-s)UW t U)\\ x \\nWt(l-Il)y( S )\\ 
Jo 

because Uy(0) = z{0) = Hz(0). Since LW tt IT < we finally have 

||r%-z)|| < C{e v + exp{-cte-")) 
and estimate ()30.b|) follows. □ 
Remark. // ITPl^n = then z' = and y is constant at leading order. 

6.3 Finite dimensional case study 

We now discuss the implications of Lemma |H1 for the different values of \ijp. We have already 
seen that a transition occurs in the definition of B e when p = /j,. Another transition happens 
when v = 0, i.e. \i — 2p. The main results of this part are summarized in Tabled 

Case when < p, < p. If < < p, we may apply Corollary El with v = and the limit 

operator B° is 

B°(n,m) = C(n,m) . 1 . l(6(n,m) ^ 0). 
7(n, m) 

In this circumstance there is no need to separate the cases when 5(n, m) ^ 0, since 

(e~M + e~ v B s ) (n, m) = e~» ( G ^ m \ + (1)) . 

\7(n,m) J 
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Hence for large values of p compared to p, the dynamics is the same for almost degenerate 
levels and for exactly degenerate levels. After a time-layer of size 0{e tl ) the system is driven 
by Eq. (|29p. namely 

3 t z = UW^Uz, z(0) = 113/(0), 
where II the orthogonal projection onto Ker A$ D KerfiJ* = Ker(\l/ e )^ om . 

Case when p = p. If p = p, we again apply Corollary El with u = p, but the form of _B° is 
slightly different, and now 8{n, m) plays a role, namely 

B £ (n,m) =B°(n,m) = C(n,m)-. -l(5(n,m) ^ 0). 

7(n, mj^ + o[n, my 

The limit equation is as in the previous case. 

Case when p < p < 2p. If p < p < 2p, we really have three different orders of magnitude in 
Eq. ([28)1 . The form for B° is now 



B°(n,m) = C{n,m) ^ n,m l l{5{n,m) + 0), 



<5(n, m) 

and the size of the time-layer is 0{e 2p ^^). The limit equation is still unchanged. 

Case when p = 2p. If p = 2p, then v = and the contribution of _B° competes with that of 
W. Therefore the projector mentioned in Corollary El is here 11^, the orthogonal projection on 
KerAjj, and the limit equation is 

d t z = U A (B° + W) t U A z, z(0) = U A y(0), 

where 



B°(n,m) = C{n,m) ]( n,m l l{5{n,m) + 0), 



S(n, m) 2 

after a time-layer of size 0{e^). In this case, and in the following one, the projection only 
depends on A (and not on B°). The constant c is simply c = min|| a; ||<i ! n A a ; =o —(A$x,x). 

Case when p > 2p or B £ = 0. If p > 2p, then v < and B e — > 0. This case is therefore 
treated in the same way as when B £ = 0. In both cases, the limit equation is 

d t z = UaW^Uaz, z{0) = U A y(0), 

after a time-layer of size 0(e^). Hence for small values of p compared to p, two almost de- 
generate levels n and m for which 5(n,m) ^ are already too far apart to resonate with the 
wave. 
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6.4 Case when fi = 0. 

In the case when p = 0, Section does not yield an interesting result since pa — * s °f order 
0(1). Therefore we have to use the transition rates obtained in Section HJ namely 

(* e )(n,m) = 2\V(n,m)\ 2 V ^ , . 7 ^' m J ^ ^|^| 2 . 

For a finite or an infinite number of energy levels, we have at leading order p^ = p^ v , where 

and 

(* >(n,m) = 2|V(n,m)| 2 £ ' |0,| 2 . 

The discussion of the long time behavior and equilibrium state for such an equation is the same 
as for Eq. (|23j) (see Appendix 18. 

6.5 Infinite dimensional case study 

In the infinite dimensional case the convergence of B £ towards B° is not sufficient to conclude. 
However, we can use series expansions when they only include one non-positive order in e, i.e. 

(^ £ ) dom + W = e- a A + o(l). 

In this case, we consider at leading order the solution p^ pp of 

d t pf p = e-^ pp pf p (31) 

where \|> ap P is homogeneous of order 0(1). In this case, time has to be changed into if = e~ u t. 
With this time scale, the equation is not singular any more. There remains to list the cases 
when such an homogeneous rate operator occurs, restricting the discussion to p > 0. As for a 
finite number of levels, the value of the ratio p/p is crucial. We once more have to consider the 
two cases p < p and p > p. The main results of this section are summarized in Tabled 

Case when p < p. If p < p, series expansions yield that 

(y £ ) dom (n,m) = C(n,m) f - . + 0(e 2p -^). 

7(n, m) 

If 2p-3p> 0, i.e. p < 2p/3 

(y £ ) dom {n,m) =: £~^ app (n, m) + o(l). 
Since v > 0, a = p and rates are homogeneous only if W = 0. 
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Case when /j = p. In the case when \i = p no series expansion is needed and we have exactly 



(* E )< i ™(n,m)=C(n,m)- „, 

7(71, mj 2 + o(n, m) 2 

Once more, we can only conclude if W = 0. 

Case when \i > p. If fi < p, series expansions lead to 



(f £ } dom (n,m) = C(n,m) 



7(71, m) 



l(5(n,m) = 0) 



+ C(7i,m) £ - 7 \: mj l(^,m)^0) + O( g 3 ^). 
o{n, my 

Two cases lead to only one term e~ a 'A + o(l), either 5(n,m) is always nonzero, /x > |p and 
= 0. In this case a = 2p — p. Or /j > 2p and a = /i. We once more need to assume that 
W = 0. 

In all the above cases, using the same type of estimates and integral formulations as in the 
proofs in Section |U and |5J we define an approximate solution p°| pp to Eq. such that 

II Pd -Pd PP !U-([0,T],/2) = 0(1), 

and \I/ app has the properties described in Section IHTT1 



h/p 



fa = 0, p > 0) 



< /i/p < 2/3 



/i/p = 1 



4/3 < fi/p<2 



li/p = 2 



a 



2p-/i 



# app (n,m) 



2|y(n, m)\ 2 J2 \'' 12 



7(n, m) 



(see caption) 



jjn, m) 



7(71, m) 



jjn, m) 2 + 5(n, m)j 



7(n, m) 
m) 2 



7(71, m) 



2 < ji/p < 00 



Table 2: Cases when the transition rates are single powers of e, for an infinite number of levels. 
The sum over /3's includes only resonant contributions in general, except when /i = where 
P E IT. 
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7 Restriction to a finite number of levels 

It is important to show that we can restrict the study to a finite number of levels since we have 
seen in Section |U] that we are able to give precise results on the time evolution of rate equations 
only in this case. Another perspective is numerical simulations, which in any case can only 
treat finite data. 

Consider a solution p to Bloch equations JIJ with initial datum p(0), and infinitely many 
quantum levels. In this section, we show how p may be approximated by p N , solution to Eq. (JTJ) 
with only a finite number N of levels. 

For this purpose, for all JV 6 N, we define n N , the projection of the space C N2 of infinite 
matrices onto the space oi N x N matrices, by 

(iT N u)(n, m) := u(n, m)l[n, m < N]. 

Then, the iV-level truncated system (^ v ) is defined from Eq. by 

e 2 d t p N (t,n,m) = —iuj £ (n,m)p N (t,n,m) + (TT N Q £ )(p N )(n,m) 



V(^,n,k] p N (t, k,m)-v(\,k,m) p N (t, n, k) 



k<N 



W) 



The initial datum p (0) is also naturally defined as 

^(0)1=^(0). 

The analysis of the previous sections shows that, in the limit e — > 0, p is approximated by 
the diagonal solution py to rate equations (|25|). Theorem also gives an approximation of p N 
by p& , which turns out to be solution to the truncated system obtained from Eq. (|25J1 : 

dt P f m = (vr^(^) dom + 7r N W) iP f {2 \ p^ 2 \0) = n N p(0). CST) 

To apply Theorem El we need Hypotheses ^ to El Under the additional condition that 
longitudinal relaxation coefficients are decaying enough at infinity (Hypothesis [Hj), we show 
that this truncation procedure is compatible with the evolution according to Bloch equations 
and rate equations: 

Lemma 10. Under Hypotheses Q] to\Q for all u,T,s> 0, there exists an integer N such that, 

if 11(1-7^(0)11,2 < v, then \\p-p N \\ L ~ {[W ) < 2v + C(e» + e 1 - 2 n, 
where C = C(T) is the constant from Theorem® and N has the form 

N = N (v,T)e~^+^ , 

with r] and N v given by Hypotheses^ and^ 

If, in addition, no resonance occurs between the wave and high energy levels (i.e. there exists 
M e N such that, when min(n, k) > M, the set {/3 G II] u(n, k) + (3 ■ uj = 0} is empty), then 
N has the form N = Nq(v,T), uniformly with respect to e. 
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Proof. Theorem|niapplies for both the infinite and the finite number of levels problems, therefore 

llp-P? ) ||^([o,TM 2 )<C(^ + 5 1 -^) 

and 

\\p N -p^ 2) h~( l0W <C(6» + 6 l ~ 2 n, 

and we only need to estimate the difference A := — p^'^ 2 \ which is solution to 

d t A = {(^ £ } dom + W) t A+ ((^ £ } dom - (%) dom ' N )tp"' {2) + (W- 7r N W) $ pf (2) . 

Thanks to the non-positiveness property of the operators associated with (\E f £ ) dom and W, an 
integral formulation leads to 

||A|| ioo([(W ) < c(||(l-7r Ar )p(0)|| i2 + || (W^-t^W^.p^IIlcoqp.t,^ 



+ || {W -n N W) iP ?®\\ L ~ {[0 ,nP), 
< C(||(l - 7T N )p(0)\\p + ||(v& £ ) dom - ^(V^ni^m^i 



+ \\W - ^WWi^nmoort ). 



The first term goes to zero as N goes to infinity simply because the initial datum p(0) is in 
I 2 . The third one reads 



sup \W(n, k)\ + sup \W(n, k)\ + sup V" \W(n, k) \ + sup V* \W(n, k)\ 

n > N k n t^N k>N n k 

< c I sup \ w ( n > k )\ + su p Yl \ w ( n > 



■ n>N k n k>N 



< CN~ K ( su P y^(l + n) x |iy(n,fc)| +^^(1 + k) K \W(n,k)\ 

< CN~ K , 



■ n>N ' k •' k>N 



thanks to Hypothesis El Thus, this term is also o(l) uniformly with respect to e as iV goes to 
infinity. 

Finally, taking into account the fact that (\l/ e ) dom is symmetric, the second term is 



2 sup £ |(^) dom (n, k)\ + sup k )\ 



■ n>N ' k " k>N 



n>N k fi; Lu(n,k)+f3-uj=0 

< Ce^N-^)^ sup J^il + n) {l+r > )N ^ J2 \M 2 \V(n, k)\< 

< (j £ -^]\f~( 1 +'n) N v _ 



n>7V fc /3; a)(n,Jfe)+/9-w=0 
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which vanishes in fact when no resonance occurs between the wave and high energy levels. Else, 
we obtain a o(l) as e goes to zero under the condition 

□ 



8 Appendix 

In this section, we give the lemmas concerning the modified relaxation operator from Eq. (J3J), 
in order to describe the dynamics of the solution. We also give a proof of the genericity of the 
small divisor estimates of Hypothesis [TJ 

8.1 Relaxation operators 

We first prove some non-positiveness properties of {W™ od )§. 

8.1.1 Continuity and non-positiveness 

Lemma 11. Let A(n, m) G l™l l m n Z£#. 

(i) Its associated operator A$ is bounded on the spaces l q , 1 < q < oo, and 

\\A iU \\ lq < \\A(n,m)\\^ii mn& Ju\\ lq . 

(ii) If in addition A(n,m) > 0, then for all positive integer N, the spectrum of the restriction 
of A$ to R N is contained in {Re A < 0} U {0}. 

(Hi) If A(n,m) > is symmetric, A$ is non-positive on I 2 , and the exponential exp(tA$) is 
well defined as an operator on I 2 when t > 0. Its norm is 1. 

Proof, (i) For q — oo or q — 1, the result is immediate. The remaining cases are obtained by 
interpolation. 

(ii) The localization of the eigenvalues of A$ is obtained via the Hadamard-Gerschgorin 
method applied to M := = A$, whose eigenvalues are the conjugates of those of Af if A is 
an eigenvalue of M, there exists an index n such that 

|A - M(n,n)\ < \M(n,m)\. 

Remarking that for m ^ n, M(n,m) = A$(m,n) > 0, and M(n,n) = —^2 m ^ n M(n,m), the 
conclusion is straightforward. 

(Hi) In the symmetric case, compute for all u G I 2 , 

(A$u,u) = A$(n, m)u{m)u{n) + A$(n, n)u{n) 2 

n mj^n n 

- ^EE M n , m )(H™)\ 2 + K™)l 2 ) - M m > n M n ? = °- 

The norm of exp(tA$) is 1 on I 2 for t > because A$ has a non-trivial kernel (see below). □ 
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In the same way as for item (i) of Lemma we have the estimate usually called Schur's 
lemma. 

Lemma 12. Let A G n V£l\ and 1 < q < oo. Then, if u G ^, m (N x N), 



A(n, k)u(k, m) — A(k, m)u(n, k) 



< 1 1 A (n, ?ri) 1 1 ;oo ; ^ nJ oo f i 1 1 ^ (ri, 1 1 ^ 



8.1.2 Asymptotic states of the rate equation 

We define the asymptotic state p associated with Eq. (jSJ) and with the initial datum p(0), the 
limit (in I 1 ), if it exists, of the solution Pd(t) to Eq. Q with the initial value p(0) as t goes 
to infinity. Such an asymptotic state is necessarily an equilibrium state (i.e. it belongs to the 
kernel of (W^™ ocl )f|)- As an example, the usual thermodynamic equilibrium p therm (corresponding 
to jy £ mod = W satisfying Eq. (JJjJ)) for the iV-level Bloch model is given by 



p therm (n) 



exp(-^) 
ELexp(-f) 



We study the kernel of operators A$ modelled on (W / £ mod )(j: the -finite or infinite- matrix 
A, written in the eigenstates basis e = (ei, e 2 , . . . ), has the property (V). 



A(m, n) = <^> A(n, m) = 0. 



(V) 



In particular, thanks to this property, a vanishing column in A^ corresponds to a vanishing line, 
and conversely. 

The kernel of A$ is linked to the spaces generated by elements of the basis e, 

El := Span{e m ; m E L}, 
when L is a finite subset of indices. 

Proposition 13. Let A(m,n) G l^l^ fl 1^1^ satisfy property (V), and A(m,n) > when 
m ^ n. In addition, suppose that there exists a decomposition of I 1 into A^-stable subspaces of 
the form El (each with finite dimension). 

Then, the restriction of A# to any such non-zero subspace with minimal size has a one- dimensional 
kernel. 

Such a minimal subset can be one dimensional, i.e. generated by a single e n . This corre- 
sponds exactly to the case when the n-th line (or column) of A$ vanishes. When the cardinality 
of L is greater or equal to two, El is characterized by: 

m < n G L <^ 3m =: m\ < • • • < m s := n, such that 

rrij G L for all j = 1, . . . , s 
and A$(rrij, ^ for all j — 1, . . . , s — 1. 
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Proof. When the dimension of El is one, the result is trivial. Now, suppose that the cardinality 
of L is at least 2. 

Denote by the restriction of A$ to El- Since the dimension of El is finite, the dimension 
of the kernel of A^ is the same as the dimension of the kernel of the transposed matrix, t Ai. 
Let u belong to this kernel. Written relatively to the basis {e n ; n G L}, the relation l A^u = 
reads 

^2 A(n,k)\u(n)= ^2 A{n,k)u{k) 

y k£L\{n} J k&L\{n} 

for all fifi. Each coefficient A(n, k) is non-negative, and thanks to the property (V), the sum 
J2keL\{n} A(n, k) is positive. Therefore, u{n) belongs to the convex hull of the other coordinates. 
Since this is valid for all n G L, all the coordinates must be equal, and the kernel of is 
generated by J] ne L e n- 1=1 

Corollary 14. Under the assumptions of Proposition for each initial datum p(0) G I 1 
satisfying Eq. (J2J), there is a unique p G I 1 HKer A$ such that for all minimal A^-stable subspace 

Since the decomposition of l l into /^-stable finite dimensional subspaces corresponds to the 
splitting of Eq. (jSJ) into a countable set of decoupled finite dimensional systems which have 
bounded solutions (thanks to the results of Lemma O on the spectrum of A$), a diagonal 
argument allows to extract a converging subsequence of pd{t)- Finally, the normalization of the 
trace on each minimal subspace ensures the uniqueness of the limit. This shows that the whole 
sequence (pd(t))t>o converges to p, which is actually an asymptotic state. 

In the purely infinite dimensional case, i.e. when no decomposition into finite dimensional 
Ajj-stable subspaces exists, the results above may break down. 

Proposition 15. Consider A(m,n) G l^l^D l^l^ satisfying property (V), A(m,n) > when 
m n, and either symmetric, or in Pauli form (relation Suppose that there exists a 

minimal A^-stable subspace E of l l generated by an infinite number of eigenstates e n . 
Then, the kernel of (the restriction to E of the operator on I 1 ) is {0}. 

Proof. Denoting Af the restriction of A$ to E, we have a one-to-one relation between elements 
of the kernel of A^ and elements of the kernel of t A^\ 

u G Ker A^ <^ l u G Ker l A^ in the symmetric case, 

^ \u(n) exp (^j^j ; n G L^j G Ker f Af in the "Pauli" case. 

Since we have the bounds 

< exp I -j^- I < exp I I < exp J 
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for all n, this correspondence preserves the summability property. 

Finally, the proof of Proposition 1131 shows that the kernel of (in I 1 ) is {0}, and this 
gives the result. □ 

Corollary 16. Under the assumptions of Proposition [731 for any initial datum p(0) with non- 
vanishing component in E, there is no equilibrium state p G Z 1 fl Ker A$ with the same trace as 
p(0) in E. 

Remark. In the symmetric case (when (\l/ £ ) dom = 0, W™ od = W), according to Lemma [771 
the I 2 -norm of the solution pd to Eq. © is decreasing in time; thus, it tends to a certain value 
r > 0. This means that pa approaches a limit cycle in I 2 belonging the intersection of the sphere 
\\Pd\\i 2 — r a> n d the hyperplane where the l l -norm is one (assuming for simplicity that there is 
no strict A^-stable subspace of I 1 ). In this case, only weak convergence (to zero) can occur. 



8.2 Diophantine estimates 

We show the genericity of Hypothesis [TJ 

Lemma 17. For all rj > and all real sequence u(n,m), there exists a constant C v > 0, such 
that for almost all value of the frequency vector uj = (u>i, . . . ,u r ), 

Va = («i, . . . , a r )eZ r \ {0}, V(n, k) G N 2 such that a ■ u + u(n, k) ^ 0, 

C 

\a ■ uj + u(n, k)\ > 



1 + \a\) r - 1 + r i(l + ny+ r i(l + k) 1 + r ' 



Proof. We follow the standard approach (see e.g. |AG91| ). Restricting u to a ball B in W, we 
show that the measure of the set of "bad frequencies" violating the inequality for all constant 
C is zero. 

For rj, c> 0, a G U \ {0} and (n, k) G N 2 fixed, set 

B Zy-= {^Bi l« • - + ^ ^ < (1 + | Q |)r-l-H, (1 + n) l + ,(l + k y + r, } ■ 

This limitates uo in the direction of a. Introducing a constant K which depends on the size of 
B only, we obtain 

Kc 

meaS i B Zn,k) < Jl 



(1 + |a|)'- 1 -H7(l + + k) 1+r > 

Now, with rj, c > fixed, the measure of the set of frequencies for which the inequality is false 
at least for some (a, n, k) is less than the sum (over a, n and k) of the ones above, and thus is 
O(c). □ 
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